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FOREWORD 


SWOE  Report  89-2,  June  1989,  was  prepared  by  Dr.  H.G.  Hughes  of  Naval  Ocean 
Systems  Center,  San  Diego,  California. 

This  report  is  a  contribution  to  the  Smart  Weapons  Operability  Enhancement 
(SWOE)  Program.  SWOE  is  a  coordinated,  Army,  Navy,  Marine  Corps,  Air  Force  and 
DARPA  program  initiated  to  enhance  performance  of  future  smart  weapon  systems 
through  an  integrated  process  of  applying  knowledge  of  the  broadest  possible  range  of 
battlefield  conditions. 

Performance  of  smart  weapons  can  vary  widely,  depending  on  the  environment  in 
which  the  systems  operate.  Temporal  and  spatial  dynamics  significantly  impact  weapon 
performance.  Testing  of  developmental  weapon  systems  has  been  limited  to  a  few  selected 
combinations  of  targets  and  environment  conditions,  primarily  because  of  the  high  costs  of 
full-scale  field  tests  and  limited  access  to  the  areas  or  events  for  which  performance  data 
are  required. 

Performance  predictions  are  needed  for  a  broad  range  of  background 
environmental  conditions  and  targets.  Meeting  this  need  takes  advantage  of  significant 
DoD  investments  by  Army,  Navy,  Marine  Corps  and  Air  Force  in  1)  basic  and  applied 
environmental  research,  data  collection,  analysis,  modeling  and  rendering  capabilities,  2) 
extensive  target  measurement  capabilities  and  geometry  models,  and  3)  currently  available 
computational  capabilities.  The  SWOE  program  takes  advantage  of  these  DoD 
investments  to  produce  an  integrated  process. 

SWOE  is  developing,  validating,  and  demonstrating  the  capability  of  this 
integrated  process  to  handle  complex  target  and  background  environment  interactions  for 
a  world-wide  range  of  battlefield  conditions.  SWOE  is  providing  the  DoD  smart  weapons 
and  autonomous  target  recognition  (ATR)  communities  with  a  validated  capability  to 
integrate  measurement,  information  base,  modeling  and  scene  rendering  techniques  for 
complex  environments.  The  result  of  a  DoD-wide  partnership,  this  effort  works  in  concert 
with  both  advanced  weapon  system  developers  and  major  weapon  system  test  and 
evaluation  programs. 

The  SWOE  program  started  in  FY89  under  Balanced  Technology  Initiative  (BTI) 
sponsorship.  Present  sponsorship  is  by  the  U.S.  Army  Corps  of  Engineers  (lead  service), 
the  individual  services,  and  the  Joint  Test  and  Evaluation  (JT&E)  program  of  the  Office  of 
the  Director  of  Defense  Research  and  Engineering  (DDR&E),  Office  of  the  Secretary  of 
Defense  (OSD). 

The  Program  Director  is  Dr.  L.E.  Link,  Technical  Director  of  the  U.S.  Army,  Cold 
Regions  Research  and  Engineering  Laboratory  (CRREL).  The  Program  Manager  is  Dr.  J.P. 
Welsh,  CRREL.  The  Integration  Manager  is  Mr.  Richard  Palmer,  CRREL.  The  task  areas 
and  their  managers  are  as  follows:  Modeling  Task  Area,  LTC  George  G.  Koenig,  USAF, 
Geophysics  Laboratory  (GL),  of  the  Air  Force  Phillips  Laboratories;  Information  Bases 
Task  Area,  Mr.  Harold  W.  West,  PE,  U.S.  Army  Engineer,  Waterways  Experiment  Station 
(WES);  Scene  Rendering  Task  Area,  Mr.  Mike  Hardaway,  Corps  of  Engineers, 
Topographic  Engineering  Center  (TEC);  Validation  Task  Area,  Dr.  Jon  Martin, 
Atmospheric  Sciences  Laboratory  (ASL)  of  the  Army  Materiel  Command. 
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INTRODUCTION 


The  spectral  radiance  of  the  sky  must  be  determined  in  order  to  predict 
the  performance  of  passive  infrared  surveillance  systems  operating  against  air¬ 
borne  targets.  Spectral  radiance  is  especially  important  when  viewing  targets 
close  to  the  horizon.  For  these  low  angles  of  viewing,  the  smooth  sea  is  a  poor 
emitter  of  infrared  wavelengths,  and  its  radiance  is  primarily  determined  by  the 
emissions  of  the  intervening  atmosphere  and  reflections  from  the  sky.  For  a 
wind-ruffled  sea,  the  emissions  from  the  individual  wave  facets  must  also  be  ac¬ 
counted  for  in  the  background  scene.  In  the  spectral  regions  about  the  6.3-^m 
water  vapor  band,  where  the  absorption,  and  consequently  the  emissivity,  is 
high,  the  sky  is  essentially  a  blackbody  at  the  temperature  of  the  lowest  atmos¬ 
pheric  layer,  with  little  variation  with  elevation  angle.  In  the  8-12-gm  region, 
the  absorption  (and  emissivity)  depends  on  the  optical  path  length,  and  therefore 
the  radiance  will  be  strongly  dependent  on  the  elevation  angle.  Cloud  type  and 
coverage  will  also  have  strong  effects  on  radiance,  depending  upon  the  elevations 
and  temperatures  of  the  clouds. 

In  the  absence  of  radiometric  measurements,  we  must  presently  rely  on 
the  LOWTRAN  6  atmospheric  transmittance/radiance  computer  code  (Kneizys,  et 
al.,  1983),  along  with  measured  meteorological  parameters  and  models  of  aerosol 
size  distributions,  to  predict  the  infrared  background  emissions.  The  LOWTRAN 
6  code  has  proven  to  be  a  versatile  tool  in  predicting  atmospheric  radiance  above 
the  horizon  (Hughes  and  Jensen,  1988).  A  model  has  recently  been  developed  of 
the  effective  radiance  of  the  sea  surface  as  a  function  of  the  viewing  zenith  angle 
and  sensor  height  (Wollenweber,  1988b),  and  has  been  incorporated  into 
LOWTRAN  6  for  processing  on  an  HP-9020  computer.  The  model  (see  Appendix 
A)  uses  the  wave  slope  statistics  of  Cox  and  Munk,  1954,  which  is  based  on  sur¬ 
face  wind  speed  and  direction  relative  to  the  look  angle.  The  model  is  currently 
limited  to  the  Navy  Maritime  Aerosol  Model  (see  Appendix  B),  which  in  addition 
to  the  standard  meteorological  parameters,  requires  as  inputs  the  current  and 
24-hour  averaged  wind  speeds  and  the  origin  of  the  air  mass.  Basically,  the 
NOSC-developed  model  calculates  the  total  contributions  to  radiance  at  the  sen¬ 
sor  from  the  emissions  of  the  intervening  atmosphere  and  the  sky  reflections  and 
surface  emissions  from  the  wave  slope  surfaces.  To-date  this  model  has  received 
very  little  validation,  and  the  computer  code  has  not  yet  been  made  available  to 
the  scientific  community. 

In  this  report,  measurements  of  near-horizon  infrared  (8-12-gm)  radiances 
made  with  a  calibrated  thermal  imaging  system  are  presented.  The  radiance 
“scenes”  (thermograms)  are  used  to  compare  the  mean  equivalent  blackbody  tem¬ 
peratures  corresponding  to  an  area  1  degree  above  and  1  degree  below  the  hori¬ 
zon  during  different  meteorological  and  surface-wind-speed  conditions.  Airborne 
measurements  of  the  vertical  profiles  of  meteorological  parameters  and  sea  sur¬ 
face  temperatures  are  then  used  to  evaluate  the  ability  of  the  radiance  algorithms 
with  the  LOWTRAN  6  code  to  predict  the  measured  results. 

MEASUREMENTS 

For  this  study,  a  Piper  Navajo  aircraft,  equipped  with  Rosemount  tem¬ 
perature  and  pressure  probes  and  an  EG&G  dewpoint  sensor,  made  vertical 
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spirals  over  the  ocean  to  obtain  temperature,  relative  humidity,  and  pressure, 
which  are  required  inputs  to  the  LOWTRAN  6  computer  code  for  calculating  the 
sea  and  sky  radiances.  A  Barnes  PRT-5  radiation  thermometer  was  also  onboard 
the  aircraft  to  measure  the  actual  sea  surface  temperatures  from  low-level 
(<20m),  constant-altitude  flights.  At  the  time  the  meteorological  parameters 

were  obtained,  measurements  of  infrared  (8-12-(im)  horizon  radiances  were  also 


780)  and  a  2.95-degree  field-of-view  lens. 


For  these  measurements,  the  scanner  was  located  at  an  elevation  of  33  m 
on  the  Point  Loma  peninsula  in  San  Diego  and  was  directed  westward  over  the 
ocean.  The  response  of  the  system  is  determined  by  placing  a  blackbody  of 
known  temperature  (±0.1°C  for  temperatures  <50°C)  in  front  of  the  lens  aper¬ 
ture.  The  digitized  video  signal  transfer  function  of  the  system  then  allows  the 
blackbody  temperature  to  be  reproduced  to  within  ±0.2  °C. 

Examples  of  measured  radiance  scenes  are  shown  in  Fig.  1  and  2  for  25 
November  and  29  September  1987.  The  data  correspond  to  low  and  moderate 
current  wind  speeds  (kc  )  of  2.5  m/s  and  7.5  m/s  for  the  respective  days.  The 
Thermal  Video  Processor  System  (THERMOTEKNIX)  available  with  the  AGA 
system  also  allows  the  thermal  scene  to  be  displayed  in  a  format  consisting  of  128 
pixel  lines  (0.023  deg/line).  The  effective  blackbody  temperature  corresponding 
to  each  pixel  can  then  be  displayed  on  the  screen  by  positioning  a  cursor  at  the 
appropriate  position.  Vertical  temperature  profiles  in  the  scene  are  displayed  on 
the  right-side  of  the  thermogram  as  shown  in  Fig.  la  and  2a.  In  each  case,  the 
horizontal  cursor  is  situated  on  the  pixel  corresponding  to  the  maximum  tem¬ 
perature  (20.4 °C  or  3.46  mW/cm*  sr  and  16.5 °C  or  3.23  mW/cm*  sr)  which  is 
taken  to  coincide  with  the  infrared  horizon.  The  increase  in  temperature  with 
altitude  (decreasing  zenith  angle)  in  the  29  September  scene  reflects  the  increase 
in  air  temperature  with  altitude  shown  in  the  Fig.  3.  Similarly,  the  decrease  in 
temperature  with  altitude  in  the  25  November  scene  can  be  related  to  the  de¬ 
crease  in  air  temperature  with  altitude.  The  larger  sky  radiance  values  on  29 
September  may  also  partially  result  from  the  lower  relative  humidity  above  200 
m  (Fig.  4),  which  decreased  the  number  of  activated  aerosols  available  to  scatter 
the  radiation. 


The  THERMOTEKNIX  system  also  allows  the  relative  frequency  of  occur¬ 
rence  of  pixels  with  a  specific  temperature  to  be  calculated  for  selected  rectangu¬ 
lar  areas  within  a  scene.  When  the  histogram  is  plotted,  a  Gaussian  curve 
(shown  with  artist-enhanced  dots)  is  calculated  which  best  fits  the  temperature 
points.  In  both  the  low  and  moderate  wind-speed  examples,  the  sky  temperature 
distributions  (Fig.  lb  and  2b)  are  more  closely  Gaussian  than  are  those  of  the  sea 
(Fig.  lc  and  2c)  which  show  a  banded  feature.  The  mean  temperature  differences 
for  the  two  examples  are  3.6°C  and  1.4 °C,  respectively.  In  the  higher  wind  speed 
(9.3  m/s)  example  shown  in  Fig.  5  for  19  March  1987,  when  whitecaps  were  pre¬ 
sent,  the  sea  temperature  distribution  (Fig.  5c)  is  remarkably  similar  to  a  Gauss¬ 
ian  distribution.  The  mean  temperature  difference  between  the  sky  and  sea  in 
this  case  is  only  0.2  °C.  In  contrast,  the  sea  and  sky  mean  temperature  distribu¬ 
tions  shown  in  Fig.  6  during  stratus  cloud  conditions  (29  April  1987)  were  found 
to  be  equal,  and  both  were  closely  Gaussian.  The  surface  wind  speed  was  not 
measured  in  this  example. 


For  this  study,  18  thermograms  recorded  during  different  meteorological 
and  wind-speed  conditions  in  1987  and  1988  were  used  to  compare  the  mean 
cloud-free  sky  and  sea  temperatures.  In  Table  1,  the  measured  mean  sky  and  sea 
temperatures  and  the  current  wind  speed  are  shown.  In  12  of  the  cases,  the 
24-hour  averaged  wind  speed  (  V)  and  the  actual  sea  surface  temperature  (tss) 
measured  by  the  aircraft  were  available.  Also  shown  are  the  combinations  of  air 
mass  factors  (AM)  (see  Appendix  B)  and  surface  visibilities  (VIS),  which  allow  the 
horizon  pixel  radiance,  calculated  using  LOWTRAN  6  (with  the  measured  pro¬ 
files  of  meteorological  parameters)  to  exactly  match  the  measured  value  (Hughes 
and  Jensen,  1988).  For  this  data  set,  the  mean  temperatures  were  highly  corre¬ 
lated  (correlation  coefficient,  r  =  0.92),  as  shown  in  Fig.  7,  with  the  sea  tempera¬ 
ture  being  less  than  the  sky  temperature,  as  indicated  by  the  dashed  line  for 
one-to-one  correspondence.  The  mean  temperatures  differed  the  most  during  low 
wind-speed  conditions.  As  shown  in  Fig.  8,  the  differences  were  found  to  de¬ 
crease  with  increasing  wind  speed. 

CALCULATION  OF  BACKGROUND  SCENES 

Earlier  work  (Hughes  and  Jensen,  1988)  has  demonstrated  that  the  sky 
radiances  measured  at  individual  pixel  lines  within  about  1  degree  of  the  horizon 
can  be  modeled  very  closely  by  using  the  LOWTRAN  6  code.  Figures  9  and  10 
show  the  comparison  of  the  measured  and  calculated  radiances  for  zenith  angles 
within  about  1  degree  below  the  horizon  for  the  scenes  in  Fig.  1  and  2.  In  the 
calculations,  the  meteorological  profiles  were  divided  into  33  layers,  as  allowed 
by  LOWTRAN  6.  The  lower  layers  of  the  profiles  were  also  divided  into  sub¬ 
layers  containing  the  same  amount  of  absorbing  and  scattering  material  and  the 
temperature  of  the  original  layer.  This  artificial  layering  has  been  found  neces¬ 
sary  (Wollenweber,  1988a)  to  remove  the  anomolous  dip  (Hughes,  1987)  which 
occurs  when  aerosols  are  included  in  the  radiance  calculations  for  zenith  angles 
close  to  90  degrees.  As  the  AGA  scanner  could  not  be  accurately  plumbed,  the 
zenith  angle  of  the  infrared  horizon  (the  pixel  corresponding  to  the  maximum  ra¬ 
diance)  in  each  thermogram  was  taken  to  be  one-half  of  a  pixel  less  than  the  an¬ 
gle  (0  =  90.179  degrees  on  29  September  and  0  =  90.174  degrees  on  25  Novem¬ 
ber)  at  which  the  LOWTRAN  calculations  indicated  the  refracted  ray  path  first 
struck  the  earth.  In  both  cases,  the  major  contributor  to  the  total  radiance  just 
below  the  horizon  is  the  path  emission  N(G)P  .  While  the  reflected  sky  radiance 
N(8)rSk  and  the  surface  emission  N(Q)SS  are  small,  their  contribution  to  the  total 
radiance  N(G)lot  at  this  low  level  (33  m)  of  observation  cannot  be  neglected.  It  is 
interesting  to  note  only  a  small  reversal  of  the  relative  magnitudes  of  the  re¬ 
flected  sky  radiances  and  surface  emissions  between  the  two  sets  of  calculations, 
which  demonstrates  the  small  influence  of  the  wave  slopes  for  the  moderate  wind 
speeds  on  29  September. 

In  a  similar  manner,  the  mean  equivalent  blackbody  temperatures  corre¬ 
sponding  to  1-degree  areas  above  and  below  the  horizon  were  calculated  for  the 
12  cases  where  vertical  profiles  of  meterological  parameters  were  available.  In 
Fig.  11,  the  measured  and  calculated  mean  sky  temperatures,  Tm(Sky)  and  Tc(Sky) , 
are  in  close  agreement  (r  =  0.99).  In  Fig.  12,  the  calculated  mean  sea  tempera¬ 
tures,  Tc(Sea),  show  slightly  more  variations  with  the  measured  values,  Tm(Sea) , 
but  the  two  are  still  fairly  well  correlated  (r  =  0.90). 
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DISCUSSION 


The  demonstrated  agreement  between  the  measured  and  calculated  mean 
temperatures  places  confidence  in  the  usefulness  of  the  models  to  calculate  radi¬ 
ances  as  viewed  from  other  altitudes.  In  a  practical  sense,  the  sea  radiance 


model  needs  to  be  improved  to  reduce  computer  running  time.  In  its  present 


quired  to  calculate  the  sea  radiance  at  one  zenith  angle  of  viewing  as  compared 
to  approximately  5  seconds  for  the  sky  radiance.  In  this  study,  approximately  36 
minutes  of  computer  time  was  required  to  calculate  the  mean  sea  temperature 
within  1  degree  of  the  horizon  on  each  day.  The  m^jor  contributor  to  the  compu¬ 
tational  time  is  the  myriad  of  sky-radiance  reflections  and  surface-emission  cal¬ 
culations  that  must  be  made  to  account  for  the  radiance  from  all  the  wave  facets. 


I«l 


Computational  time  could  be  saved  if  the  surface  emissivity  and  reflectivity  cal¬ 
culations  were  made  at  a  single  representative  wavelength  instead  of  having  to 
be  averaged  over  a  wavelength  band  for  each  wave  slope.  Another  way  of  reduc¬ 
ing  the  computational  time  is  to  develop  an  empirical  relationship  to  directly  re¬ 
late  the  apparent  mean  sea  temperature  to  the  calculated  mean  sky  temperature. 


Inherent  in  the  observed  variations  between  the  mean  sky  and  sea  tem¬ 
perature  differences  in  Fig.  5  are  the  surface  wind  speeds  and  actual  sea  surface 
temperatures.  If  these  parameters  are  included  in  a  multiple  regression  analysis 
for  three  independent  variables,  we  obtain  the  following  relationship: 


Tm(Sea)  =  1.09  Tm(Sky)  +  0.37  Vc  +  0.24  Tss  -  10.46  (1) 


The  correlation  coefficient  for  the  multiple  regression  analysis  is  0.97.  Only  12 
samples  were  used  in  the  analysis  (as  compared  to  18  samples  in  the  linear  re¬ 
gression  analysis  shown  in  Fig.  5),  which  indicates  a  definite  improvement  in  the 
correlation  when  the  actual  sea  surface  temperatures  and  wind  speeds  are  in¬ 
cluded. 


Based  on  earlier  work  (Hughes  and  Jensen,  1988),  which  demonstrated 
that  the  sky  radiances  measured  at  individual  pixel  lines  within  about  1  degree  of 
the  horizon  can  be  modeled  very  closely  by  using  the  LOWTRAN  6  code,  and  the 
excellent  agreement  between  the  calculated  and  measured  mean  sky  radiances  in 
Fig.  11,  the  replacement  of  Tm(Sky)  by  a  calculated  value  of  the  mean  sky  tem¬ 
perature,  Tc(Sky) ,  would  appear  justified.  This  would  allow  the  mean  sea  tem¬ 
perature  to  be  directly  determined  by  using  Eq.  1  and  measured  profiles  of 
meteorological  parameters.  Additional  data  are  needed  to  verify  that  such  a  rela¬ 
tionship  for  predicting  the  apparent  mean  temperature  of  the  sea  is  meaningful 
for  other  meteorological  conditions,  sea  states,  and  altitudes  of  observation. 

Such  a  relationship  would  be  a  valuable  inclusion  into  a  system  designer’s  hand¬ 
book  as  a  tool  for  predicting  infrared  sea  backgrounds  in  locales  where  only  mete¬ 
orological  data  are  available. 


Appendix  A 

MATHEMATICAL  FORMULATION  OF  SEA  RADIANCE 

MODEL 

Consider  the  atmosphere  to  be  composed  of  a  number,  n,  of  isothermal 
layers  characterized  by  temperature  t)  and  transmittance  along  the  opti¬ 

cal  path  traversing  the  ith  layer  at  angle  ji,  and  v  is  the  spectral  wave  number. 
From  KirchofPs  law,  the  radiance  of  the  ith  layer  is 

N(v,i,ii)sk  =  [1  -  ra(v,i,n)]W(T,)/ir  (A-l) 

where  ra  (v,i,ji)  is  the  absorption  transmittance  and  W(Ti)  is  Planck’s  blackbody 
radiation  formula.  Then  the  spectral  radiance  reaching  the  sea  surface  through 
the  intervening  atmosphere  is 

=  [1-  r.(v,i,n)]^n  t(v,j,^W(Ti)h  (A-2) 

Summing  the  contribution  from  all  layers,  the  spectral  radiance  at  the  sea  sur¬ 
face  is  then 


N(v,nU=  S{[1-  r4vU))  nr(v,j,n)  W(Ti)h} 

1=1  L/  =  1 


(A-3) 


The  radiance  is  allowed  to  strike  a  wave  facet  on  the  ocean  surface  with  a  Gauss¬ 
ian  distribution  (Cox  and  Munk,  1954)  of  angular  tilts  a  and  p  in  the  up-wind 
and  cross-wind  directions,  respectively,  so  that  an  amount  is  reflected 

into  the  sensor  at  an  altitude  within  the  mth  layer.  The  probability  that  radi¬ 
ance  hits  the  facet  is  equal  to  the  probability  that  the  wave  slope  exists,  i.e., 

Niv^Y  =  P(Sx,Sy)NMsk  (A-4) 

where 

P(Sx,Sy )  =  l/(2o^)EXP[0.05(Si/o5  +  S^)]  (A-5) 

and  Sx  =  tana,  =  tanp,  =  0.003  +  1. 92x10s VCy  ciy  =  3. 16x103 Vc,  with  Vc 
being  equal  to  the  current  wind  speed  in  the  azimuthal  direction  <J>  with  respect 
to  the  sensor.  Then  the  total  spectral  radiance  that  is  reflected  from  all  the 
wave  facets  into  the  line-of-sight  of  the  detector  located  in  the  mth  layer  is 

[m-  1  "1  2  R 

n  r(j,d)  p  (v,Si)P(Sx  ,  Sy)N(v,ik)rsk  (A-6) 

where  R(v,Ci)  is  the  complex  reflectivity  of  seawater  at  the  reflection  angle  fl.  In 
the  above  equations,  both  and  n  are  implicit  functions  of  Sx  and  Sy  given  by 
(Wollen  weber,  1988b) 
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cosp.  =  (2S*44)cos0'cos4>  +  (25^,44)0080' cos4>  -  (B44)sin0' 
cosfl  =  (5x/A)cos0'cos<j))  +  (5^44)0080' cos<J>  +  (144)sin0' 


(A-7) 


(A— 8) 

where  A  =  S*  +  Sy  +  1,  B  =  Si  +  - 1,  and  6'  is  the  sensors’  zenith  angle  at 

the  sea  surface  reflection  point, 


In  a  similar  manner,  the  spectral  radiances  emitted  by  the  sea  surface 
wave  facets  (Nss)  and  the  path  radiance  (/Vp)  which  reach  the  sensor  at  the  zenith 
angle  0  are  given  by 


N(v,Q)ss  = 


*£(SX  ,5y)[l-i?(v,n)]W(TM)/TT 


(A-9) 


and 


m 

N(v,Q)P  =  X[l-ra(i),t, 

I  =  1 

where  Tss  is  the  sea  surface  temperature,  and  again,  the  angle  n  is  implicit  in  the 
reflection  angle  fl. 

Then,  the  totcil  spectral  radiance  reaching  the  detector  is  the  sum  of  the 
three  components: 

N(v,Q)tot  =  N(v,Q)rsk  +  N(v,Q)ss  +  N(v,Q)p  (A-ll) 

The  total  spectral  radiance  must  then  be  averaged  over  the  response  of  the  sensor 
system,  which  in  this  case  is  taken  to  be  the  8-  to  12-^m  wavelength  band.  Sub¬ 
routines  have  been  introduced  into  LOWTRAN  6  (Wollenweber,  1988b)  to  calcu¬ 
late  the  band-averaged  radiance.  The  reflection  and  zenith  angles  are  calculated 
with  Eq.  A-5,  A-7,  and  A-8,  corresponding  to  the  incremented  values  of  wave 
slopes  in  the  intervals  -3oxy  <  Sxy  <  SoXJ . 


]  n  t(vJ,6) 


W(Tt)h 


(A- 10) 


Appendix  B 

DESCRIPTION  OF  THE  LOWTRAN  6  NAVY  MARITIME 

AEROSOL  MODEL 

The  particle  size  distribution  model  (at  radius  r)  is  the  sum  of  three  log¬ 
normal  distributions  given  by 
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n(r)  =  X  A  i  exp[-(ln  F]  (cm-3  •  jim-l) 

i  =  l  Jri 

(B-l) 

where 

A,  =  2000(AMF 

(B-2) 

A2  =  5.866(F  -  2.2) 

(B-3) 

A3  =  10  (0.06VC  -  2.8) 

(B-4) 

Component  A,  represents  the  contribution  by  continental  aerosols.  AM  is 
an  air  mass  parameter  that  is  allowed  to  range  between  integer  values  of  1  for 
open  ocean  and  10  for  coastal  areas  and  is  given  by 

AM  =  INT(Rn/A)  +  1  (B-5) 

where  Rn  is  the  measured  atmospheric  radon  content  expressed  in  pCi/m3.  In 
the  absence  of  radon  measurements,  the  air  mass  factor  can  be  related  to  the 
elapsed  time,  T(days),  for  the  air  mass  to  reach  the  point  of  observation: 

AM  =  J/VT[9exp(-T I  A)]  +  1  (B-6) 

Components  A2  and  A3  represent  equilibrium  sea  spray  particles  generated 
by  the  surface  wind  speed  averaged  over  24  hours  (V,  in  m/s)  and  the  current  sur¬ 
face  wind  speed  (Vc,  in  m/s),  respectively.  (It  should  be  noted  that  the  current 
wind  speed  component  is  different  from  the  value  published  in  LOWTRAN  6. 

This  modification  was  found  to  be  necessary  in  order  to  match  previously  pub¬ 
lished  measurements  of  infrared  sky  radiances  and  near-surface  aerosol  size  dis¬ 
tributions  (Hughes,  1987)).  In  Eq.  B-l,  r, ,  the  modal  radius  for  each  component 
referenced  to  a  relative  humidity  of  80%  (r,  =  0.03  p.m,  r2  =  0.24  jim,  and  r3  = 
2.0  p.m),  is  allowed  to  grow  with  relative  humidity  (RH)  according  to  the  formula 
(Fitzgerald,  1975) 

f=  K2  -  i?H/100)/6(l  -  RH/  100)1V3  (B-7) 

The  contribution  to  the  total  extinction  or  absorption  by  each  component 
can  then  be  written  as 

oe,a  (k)i  =  ( SF){ct  /  (^r.w)expr-(ln -Lyirklr}  (B-8) 

r  Jr‘ 

where  C,-  =  (O.OOIW/M,. 
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The  factor  fx  in  the  expression  for  C,  ensures  a  constant  total  number  of 
particles  as  the  relative  humidity  increases.  Qe.a  ( \,r,m )  is  the  cross-section  for 
either  the  extinction  or  absorption  normalized  to  the  geometrical  cross-section  of 
the  spherical  particle,  and  m  is  the  complex  refractive  index,  which  is  allowed  to 


change  from  that  of  dry  sea  salt  as  the  particle  deliquesces  with  increasing  hu- 


rameter  ot.a  (h)/  C,  at  discrete  wavelengths  for  four  relative  humidities  (50%, 
85%,  90%,  and  99%),  from  which  the  average  extinction  for  a  specific  wavelength 
band  and  relative  humidity  can  be  readily  determined  by  interpolation.  When  an 
observed  surface  visibility  ( VIS0 )  is  available  as  an  input  to  the  model,  the  ampli¬ 
tudes  of  the  three  components  are  adjusted  by  a  scaling  factor  ( SF)  so  that  the 
calculated  aerosol  extinction  coefficient,  °c,  at  a  wavelength  of  0.55  pm  is  the 
same  as  the  observed  extinction,  o0,  determined  from  the  relationship 


VIS0 


3.912 
a0  +  or 


where  ar  is  the  Rayleigh  contribution  to  extinction  at  0.55mm. 


(B-9) 
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Thermograms  of  near-horizon  infrared  radiances  on  25  November  1987  during  low 
wind  speed  conditions  (2.5  m/s),  (a)  Vertical  profile  of  effective  biackbody  tempera¬ 
ture  near  the  horizon,  and  temperature  histograms  corresponding  to  areas  1  degree 
above  (b)  and  1  degree  below  (c)  the  horizon.  In  each  histogram  a  Gaussian  curve 


Is  displayed  as  a  series  af  dots  which  best  lit  tbe  temperature  points. 


Thermograms  of  near-horizon  infrared  radiances  on  29  September  1987  during  mod¬ 
erate  wind  speed  conditions  (7.5  m/s)  without  the  presence  of  whitecaps.  (a)  Verti¬ 
cal  profile  of  effective  biackbody  temperature  near  the  horizon,  and  temperature 
histograms  corresponding  to  areas  1  degree  above  (b)  and  1  degree  below  (c)  the 
horizon.  In  each  histogram  a  Gaussian  curve  is  displayed  as  a  series  of  dots  which 
best  fit  the  temperature  points. 

Profiles  of  air  temperature  measured  with  altitude  on  29  September  and  25  Novem¬ 
ber  1987  off  the  coast  of  San  Diego,  California. 

Profiles  of  relative  humidity  measured  with  altitude  on  29  September  and  25  Novem¬ 
ber  1987  off  the  coast  of  San  Diego,  California. 

Thermograms  of  near-horizon  infrared  radiances  on  19  March  1987  during  moderate 
wind  speed  conditions  (9.2  m/s)  when  whitecaps  were  present,  (a)  Vertical  profile  of 
effective  biackbody  temperature  near  the  horizon,  and  temperature  histograms  corre¬ 
sponding  to  areas  1  degree  above  (b)  and  1  degree  below  (c)  the  horizon.  In  each 
histogram  a  Gaussian  curve  is  displayed  as  a  series  of  dots  which  best  fit  the  tem¬ 
perature  points. 

Thermograms  of  near-horizon  infrared  radiances  on  29  April  1987  during  stratus  cloud 
conditions  (wind  speed  not  measured),  (a)  Vertical  profile  of  effective  biackbody 
temperature  near  the  horizon,  and  temperature  histograms  corresponding  to  areas  1 
degree  above  (b)  and  1  degree  below  (c)  the  horizon.  In  each  histogram  a  Gaussian 
curve  is  displayed  as  a  series  of  dots  which  best  fit  the  temperature  points. 

Regression  analysis  of  measured  mean  sky  temperatures,  Tm(Sky)<  and  mean  sea 
temperatures,  Tm(Sea >  The  dashed  line  corresponds  to  a  one-to-one  correspon¬ 
dence  between  the  two  parameters. 

Regression  analysis  of  the  difference  between  the  measured  mean  sky  temperatures, 
Tm(Sky),  and  mean  sea  temperatures,  Tm(Sea >  and  the  surface  wind  speed. 

Comparison  of  the  measured  and  calculated  IR  radiances  for  zenith  angles  about  1 
degree  below  the  horizon  on  25  November  1987. 

Comparison  of  the  measured  and  calculated  IR  radiances  for  zenith  angles  about  1 
degree  below  the  horizon  on  29  September  1987. 

Correlation  between  the  measured,  Tm(Sky).  and  calculated,  Tc(Sky) .  mean  sky 
temperatures. 

Correlation  between  the  measured,  Tm(Sea) .  and  calculated,  Tc(Sea)  .  mean  sea 
temperatures. 
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TABLE  1.  Mean  sky,  Tm(Sky) ,  and  sea,  Tm(Sea),  temperatures,  current  wind  speed,  Vc, 

24-hour  averaged  wind  speed,  V,  and  the  actual  sea  surface  temperature,  Tss,  measured  on 
the  days  indicated.  VIS  and  AM  are  the  inferred  visibility  and  air  mass  factor,  respectively 
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